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Abstract 
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growth) and growth of slits. We show that the Lowner differential equation for a 
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1 Introduction 



Growth problems of Laplacian type (such as Hele-Shaw viscous flows) refer to dynamics 
of a moving front (an interface) between two distinct phases driven by a harmonic scalar 
field. These essentially nonlinear and non-local problems attract much attention for quite 
a long time [H [2] . 

Remarkably, Laplacian growth with vanishing surface tension possesses an integrable 
structure uncovered in [3] . Since evolution of planar domains is most naturally described 
by time-dependent conformal maps, there is no surprise that this structure is actually im- 
manent for general conformal maps and classical boundary value problems. Specifically, 
in [1] it has been shown that evolution of conformal maps is governed by an integrable 
hierarchy of nonlinear partial differential equations which is a zero dispersion version |5] 
of the 2D Toda hierarchy [6] . In fact the Lax equations for this hierarchy can be derived 
from the classical theory of conformal maps depending on parameters. 

In the present paper we have tried to give a short review of these and related devel- 
opments in the context of the integrable systems theory. Taking the Lax representation 
of the dispersionless 2D Toda (dToda) hierarchy as a starting point, we follow how it 
induces the contour dynamics and dynamics of conformal maps. The Lax function, which 
is supposed to be univalent in some neighborhood of infinity, is interpreted as a conformal 
map from a fixed reference domain to the complement of the growing domain. 

Depending on the type of solutions to the dToda hierarchy, the growing domains can 
be either smooth or singular like cuts or slits. Hence the system of dispersionless Lax 
equations serves as a master dynamical equation not only for Laplacian growth but also 
for growth of slits. Conformal maps of parametric families of slit domains are known 
to satisfy a differential equation proposed by K.Lowner in 1923 Its connection with 
nonlinear integrable equations was pointed out in |8], see also [9l [TOl [11]. Following 
|12] . we demonstrate how the radial Lowner equation arises in the context of the dToda 
hierarchy. 

The unified treatment of the two types of growth processes in the framework of 
the Toda integrable system, which we emphasize in this paper, seems to be especially 
promising in the light of the stochastic Lowner evolution (SLE) approach [T3|. This 
might give a hint how to incorporate, in an intelligent way, a stochastic ingredient into 
growth problems of Laplacian type. 

An instructive representation of solutions to the dToda hierarchy is provided by the 
large N limit of certain matrix integrals or their eigenvalue versions (the Dyson gas 
representation), with the associated growth processes being simulated by evolution of 
support of eigenvalues. In the last section we outline the Dyson gas representation for 
Laplacian growth and growth of slits. 

Section 2 contains the necessary material on the dToda hierarchy and its Lax rep- 
resentation. In section 3 we present the general solution to the hierarchy in terms of 
canonical transformations and distinguish the classes of non-degenerate and degenerate 
solutions. In section 4 we associate a contour dynamics with any solution to the Lax 
equations. A particular subclass of degenerate solutions is studied in section 5, where 
the Lowner equation is derived from the Lax equations. Finally, section 6 contains the 
Dyson gas representation for both non-degenerate and degenerate solutions. 
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2 Lax representation for the dToda hierarchy 



This section contains a standard material which we present in a form convenient for 
our purposes. For a more complete account of dispersionless hierarchies, their algebraic 
structure, solutions and applications see [5].[T^-p!7]. 

Dispersionless Lax equations. We start with the Lax representation of the dToda 
hierarchy with certain reality conditions imposed. The main object is the Lax function 
z{w) represented as a Laurent series of the form 

z[wj =rw + aQ-\ \ + ... 

w 

The leading coefficient r is assumed to be real while all other coefficients Oj are in general 
complex numbers. All the coefficients depend on deformation parameters (or "times" ) to 
(a real number) and ti, t2, ts, . . . (complex numbers) in accordance with the Lax equations 

^ = {A,{w), z{w)} , ^ = -{A,iw-'), z{w)} (1) 

where for any two functions of w, to we set 

{/ g} ■= - (2) 

' d log w dto dto d log w 



Here and below the bar is complex conjugation and f{w) means f{w). The reality 
condition thus implies that the second half of the Lax equations (with t^-derivatives) 
is obtained from the first one by complex conjugation with w on the unit circle. The 
generators of the flows are constructed as follows: 

Ak{w) = (z^iw))^ , Ao{w) = logw 
For the dToda hierarchy, the (. . .)4_-operation is 

Hereafter, {■ ■ ■)s nieans taking the terms of the Laurent series with degrees belonging to 
the subset 5* G Z (in particular, (. . .)o is the free term). Note that at k = equations 
([1]) become tautological identities. The second Lax function of the dToda hierarchy is 
z{w~^). The reality conditions (i.e. the requirement that its coefficients are complex 
conjugate to those of the z{w)) imply that it obeys the same Lax equations. It should 
be noted that is just a compact form of writing evolution equations with respect to 
real deformation parameters = TZetj^, t\ = Imtk- 

= {Ak{w) - Ak{w~^), z{w)] , = I [Ak{w) + Ak{w-^), z{w)] 

We are especially interested in the class of solutions such that z{w), for all tk in 
an open set of the space of parameters, is a univalent function in a neighborhood of 
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infinity including the exterior of tlie unit circle. Tliis means tliat in tliis neigfiborliood 
z{wi) = z{w2) if and only if Wi = W2- From now on, we assume that z{w) belongs to 
this class. In this case z{w) is a conformal map from the exterior of the unit circle to a 
domain in the complex plane containing infinity while z{w~^) is a conformal map from 
the interior of the unit circle to the complex conjugate domain. 

Let w{z) be the inverse function to the Lax function z{w). In terms of the inverse 
function, the evolution equations ([1]) acquire a simpler form: 

d log w{z) _ dAk d log w{z) _ _dAk 
dtk dto ' dik dto 

Here = Ak{w{z)), Ak = Ak{l/w{z)) are regarded as functions of z, and the derivatives 
are taken at constant z. 

By purely algebraic manipulations, one can show [5] that the compatibility conditions 
for the Lax equations ([1]) read 

dt^Akiw) - dt,A,{w) + {Ak{w), A,{w)} = 

dt^Akiw-') + dt,A,{w) + {Akiw-'),Aj{w)} = 

which is a dispersionless version of the "zero curvature" representation. Treating A^s as 
functions of z, one can rewrite them in the form similar to ([3]): 



dAj _ dAk dAj _ dAk 



dtk dtj dtk dtj 
Note that at j = this system coincides with ([3]). 



(4) 



Generating form of the Lax equations. Using a generating function of the poly- 
nomials Ak, the infinite hierarchy ([T]) can be represented as a couple of "generating 
equations" . The generating function is defined as 



-k 



7 - iz^iw)] 



log 1-^'"' 



z{wi) 



where wi = w{zi). To separate the polynomial part and the free term, we write 

/ z(w) \ , \ , rw(zi) , z(wi) — z(w) 
log 1 - = log 1 - + log^^ + log^j^ Y- 

and notice that the expansion of the first (third) term contains only positive (respectively, 
negative) powers of w while the rest is just the free term. Therefore, 

E '-^Mw) = - log(«^(^i) -w) + l log (5) 
^1 ^ 2 r 



Similarly, 



E ^^Mw~') = - log(«^(^i) - w-') + I log (6) 
^1 ^ 2 r 
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and the generating Lax equations read 



D{zi)z{w) = - I \og{w{zi) - w) + J log — , z{w) \ (7) 
I 2 w{zi) J 

- f 1 r 1 

D(zi)2;(w) = < log(w;(^i) - w"^) + - log 



2 ""wizi)' ' ' j 
where we have introduced the differential operators 

k>i k>i 

Expanding both sides in powers of Zi, one recovers eqs. ([T]). In terms of the inverse 
function, the generating Lax equations acquire a more transparent form (cf. ([3])): 

D{zi)\ogw{z2) = -dialog {w{zi) -w{z2)) + l-dt^log^^^^^ (9) 

D{z,) \ogw{z2) = dt, log (w{z,) - w-\z2)) - \dt, log ^ (10) 
Tending 2:2 — > cxd, we obtain the useful equations 

D{z) logr = -]^dto log {rw{z)) , (11) 

L'(z) logr = -^dto log (rw(^)) (12) 

Plugging them back into eqs. (Q, ( ITOi) . one can represent the latter relations in a slightly 
more compact form: 

D{zi) \og{rw{z2j) = -dtg log {rw{zi) - rw{z2)) (13) 
D{z,) \og{rw{z2)) = dt, log f 1 - I ] (14) 

\ w{Zi)w{Z2) J 

The construction of A^'s implies that the expansion of Ai:{w{z)) in a Laurent series in 
z looks like A^. = z'^ + 0(1). Furthermore, the compatibility conditions (jl]) and relation 
allow one to represent the coefficients in the form 

Ao{w{z)) = \ogw{z) = -]-dt,vo -Y^^z-' (15) 

^ k>i 

A,{w{z)) = z^-\d,^v,-Y.^-^z-\ j>l (16) 
^ k>i ^ 

A,{w-\z)) = \d,^v, + Y.^-^^-'^ ^■>1 (17) 
^ fc>i 

where are functions of the times such that dt-Vk = dt^Vj, dt Vk = dt,,Vj. The latter 
conditions allow one to introduce a real- valued function JF via = dtf.J-'. Equations 
f lT^ . ffT^ then become the dispersionless Hirota equations for the function JF. 
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3 General solution to the Lax equations 



A general solution to the differential equations ([H) is available in an implicit form [5] . To 
present it, we need an extended version of the Lax formalism. 

The idea is as follows. By the definition of the Poisson bracket, logw and to form a 
canonical pair: {logw, to} = 1- The evolution according to the Lax equations can be 
regarded as a t^-dependent canonical transformation from the pair (log if, to) to another 
canonical pair whose first member is log z{w). It is quite natural to introduce the second 
member which we denote by M. Depending on the situation, we shall treat it either 
as a function of z and to or as a function of w and to through the composition M = 
M{z{w,to),to) (it also depends on the deformation parameters t^). To find what is M, 
we note that the condition {log 2;, M} = 1 can be identically rewritten as dtgM{z) = 
zdz log w{z, to). This determines M up to a term depending only on z. The latter is fixed 
if one requires M to obey the same Lax equations ([I]). To wit, equation dt,^M = {Ak, M} 
(where the derivatives are taken at constant w) is equivalent to 

dt,M{z) = wd^Akdt,M{z) = zd.Ak 

Taking into account ( |T5|) . ( |T6i) . we can write 

M = J2 khz\w) + to + ^ VkZ-\w) (18) 

fc>l k>l 

It is the quasiclassical (dispersionless) limit of the Orlov-Shulman operator [18]. Its 
geometric meaning depends on the choice of a particular solution. In a similar way, one 
can construct the conjugate Orlov-Shulman function, M{z), such that the transformation 
{logw, to) {log z~^ {w~^), M{z{w~^)) is canonical and M obeys the same Lax equations. 
The Lax equations imply that the composition of the canonical transformations 

(log2;,M) (log I/;, to) ^ {log z'\M) 

does not depend on t^, i.e., it is an integral of motion. Moreover, any t^-independent 
canonical transformation {log z,M) —>■ {log z~^,M) between the Laurent series of the 
form prescribed above generates a solution to the dToda hierarchy. A detailed proof can 
be found in [5]. 

More precisely, let {log f{w,tQ), g{w, to)) be a canonical pair: {log f,g} = 1. Sup- 
pose that the functions z, z, M, M of the form as above are connected by the functional 
relations 

l/z{w-') = f {z{w), M{z{w))) , M{z{w-')) = g {z{w), M{z{w))) (19) 

Then the function z{w) obeys the hierarchy of the Lax equations and its coefficients (as 
functions of t^'s) thus obey the dToda hierarchy. Conversely, any solution of the dToda 
hierarchy admits a representation of this form with some (/, (yf)-pair. Note that the reality 
conditions imply the following constraints on the functions /, g: 

{r\w,to), g{w,to)) =w, g[f-\w,to), g{w,to)) =to (20) 

which is a sort of the "unitarity condition" for the canonical transformation (logw, to) 
{log f{w, to), g{w, to)). 
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This construction can be made more explicit by introducing the generating function 
of the canonical transformation (log w, to) {^og f,g). An important class of solutions 
corresponds to the canonical transformations (log z, M) (log^~^, M) defined by means 
of a generating function U{z, z) [19] : 

M = zd,U{z, z), M = zdsUiz, z) (21) 

Here t/(2, z) can be an arbitrary differentiable real- valued function of z, z. This form of 
the canonical transformation implies that the functions z[w) and z{w~^') are algebraically 
independent. (By algebraic dependence we mean here existence of a t^-independent 
function of two variables R{z,z) such that R{z{w), z{w~^)) = for all tk-) These are 
solutions of generic type. We call them non-degenerate. For non-degenerate solutions 
the "string equation" 

where Uzz{z,z) = dzdzU{z,z) holds true. It is obtained by plugging M from fl2T]) into 
the canonical relation {z, M} = z. 

The origin of the string equation can be understood in a simpler way as follows. 
From the Lax equations ([I]) and the Jacobi identity for the Poisson bracket it follows 
that {z{w),z{w~^)} obeys the same Lax equations: 

dt,{z{^)rz{w-^)} = {Akiw), {ziw),ziw-')}} 

Therefore, any relation of the form 

{z{w), z{w-^)} = uj{z{w), z{w-^)) (23) 

where uj{z, z) is an arbitrary tfc-independent function of two variables such that io{z, z) = 
uj{z, z) is consistent with the hierarchy. The approach based on the canonical transfor- 
mations clarifies the meaning of this function and makes it clear that in fact any solution 
obeys a string equation of the form fl23|) . 

Canonical transformations such that the functions z{w) and z{w~^) appear to be 
algebraically dependent can not be represented in the form (pTI) . They correspond to 
solutions which we call degenerate. For degenerate solutions, the Poisson bracket {z, z} 
vanishes. Conversely, the relation {z, z} = imphes the algebraic dependence. Indeed, 
in terms of the function z{w~^{z)) the Lax equation for z{w~^) reads: 

dt^z{w^^{z)) = dzAk {z, z} 

(the derivatives are taken at constant z), so at {z,z} = we have dti^z{w^^{z)) = for 
all tfc- This just means that the relation between z{w) and z{w~^) is t^-independent. A 
particular subclass of degenerate solutions, together with their geometric interpretation, 
is discussed below in section 5. 

4 Contour dynamics 

The Lax equations ([T]) can be understood as equations of a contour dynamics. The 
contour is the image of the unit circle, i.e., z{e^^), < 6 < 2it. Let us call it the Lax 
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z(w) 




\w\=l 

Figure 1: The Lax contour. 

contour (Fig. [T]) and denote it by 7. It depends on the deformation parameters according 
to the Lax equations. 

We need a general kinematic relation. Let {x{a,t),y{a,t)) be any parameterizations 
of a moving contour in the plane, then the normal velocity of the contour points is 

K = ^ {d^xdtv - d^ydtx) 



where dl = y {dxy + {dyy is the line element along the contour. 

Applying this formula to the Lax contour z(e*^) with the specific parametrization 
a = 6 and t = to with all other t^s fixed, we get the normal velocity of the Lax contour 
7 at the points z{w), \w\ = 1: 

_ {z{w), z{w-')} 

Here z'{w) = dwz{w) and the Poisson bracket in the numerator is given by ([2]). More 
generally, the normal velocities corresponding to the changes of the (real) times t^ = 
TZetk, t], = Xmtk are given by 

t'W|''' (^--W-^''''"''>) 

where 

(pk{w) ■= wd^Ak{w) (27) 

The function z{w) provides a time- dependent conformal map from the exterior of the 
unit circle onto the exterior of the Lax contour. 



Non-degenerate solutions. The non-degenerate solutions corresponding to the canon- 
ical transformation with the generating function U{z, z) fl^Tl) have a clear interpretation 
in terms of contour dynamics. Eq. together with the string equation (!22|) states 

that the normal velocity of the Lax contour at the point 2; G 7 is equal to 
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Eqs. f|T8|) . fl2Tl) allow us to express the deformation parameters in terms of the moving 
contour: 



z"^-\w)M{z{w))dz{w) = ^— (f z-''d,Udz, k>l (29) 



27lik j'|u;|=l 2TTik 

We stress that ti, ^2, • • • are kept constant, so they are integrals of motion for the contour 
dynamics fl28|) . 

A particularly important case is U {z, z) = zz which corresponds to the canonical 
transformation z = z~^M, M = M (i.e., M = M = zz). In this case the normal velocity 
is given by 

K(^) = -^K(^)|, ze^ (31) 

Note that |w'(-2)| is equal to the normal derivative dn\og\w{z)\ of the solution to the 
Laplace equation with a source at infinity and the Dirichlet boundary condition on the 
contour. We thus see that (|3T|) is identical to the Darcy law for the dynamics of interface 
between viscous and non-viscous fluids confined in the radial Hele-Shaw cell, assuming 
vanishing surface tension on the interface. Formulas ( 129|) state that tk = 2^ ^7 z~^zdz 
are harmonic moments of the exterior of the contour 7. Their conservation in the course 
of the Laplacian growth dynamics was first established by S.Richardson [20]. Eq. (!29l) 
states that the time variable to should be identified with area (divided by vr) of the interior 
domain encircled by 7. 

The Laplacian growth with the source at a finite point zq corresponds to the same 
function f/(z, z) = zz and the vector field dt^ + D{zq) + D{zq) in the space of deformation 
parameters. Indeed, using eqs. fl2S]) . fl2B]) and relations ([S]), ([S]) we obtain the normal 
velocity 

, . dnG{z, Zq) , , 

where 

G{z, zo) = log 



1 - w{z)w{zo) 

is the Green function of the Dirichlet boundary value problem in the exterior of the Lax 



contour. 



Degenerate solutions. Degenerate solutions describe evolution of singular contours 
like growth of slits or cuts in the plane. Since the Poisson bracket {z, z} vanishes, it 
might seem from that the velocity of the contour vanishes as well and so there is no 
dynamics at all. In fact this is not exactly the case: Vn does vanish unless z'{w) = 0. We 
see that the growth is possible only at the points that are images of the critical points 
of the conformal map z{w) lying on the boundary (on the unit circle). This means that 
only endpoints of arcs can move while other boundary points remain fixed. We see that 
the Lax contour or at least a finite part of it degenerates into an arc of a fixed curve F 
swept twice (back and forth), or into a collection of such arcs, and the evolution consists 
in moving the endpoints of the arc along the same curve F (see an example in Fig. [2]). 
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Figure 2: A degenerate Lax contour. The arc can move along a fixed curve T. 



This agrees with the fact that for degenerate solutions the functions z{w) and z{w~^) 
are algebraically dependent: the t^-independent relation R{z{w), z{w~^)) = is just the 
equation of the curve F. The deformation parameters do not admit so transparent 
interpretation as in the non-degenerate case. 



5 Reductions of rank 1 and radial Lowner equation 

Any particular solution of the dToda hierarchy can be regarded as a transition from the 
variables to,ti,t2, ■ ■ ■ to the variables r, ao, ai, . . . which are coefficients of the Lax function 
z{w): z{w) = rw + aQ + aiZ"^ + . . .. The non-degenerate solutions result in true changes of 
variables, i.e., the Jacobian of this transition does not vanish. For degenerate solutions, 
the Jacobi matrix is degenerate. In this section we present a detailed analysis of the 
simplest nontrivial case when this matrix is of rank 1. We say that the corresponding 
solutions are reductions of rank 1 (of the dToda hierarchy). As is easy to see, the reduction 
of rank 1 implies that the Lax function and thus w{z) depends on all the times tj through 
only one independent function q = q{{tj}): 

z{w; {tj}) = z{w, q) , w{z] {tj}) = w{z, q) (33) 

Our goal is to characterize the form of the function w{z) = w{z,q) or z{w) = z{w,q) 
consistent with the dToda hierarchy. Without loss of generality, we set q := logr. We 
shall see that the consistency condition is the radial Lowner equation [7] 

dwjz) r]{q)+w{z) 

-j; — = TT (34) 

oq ri[q) — w[z) 

where r]{q) is arbitrary continuous function of q such that \r]{q)\ = 1 (the "driving func- 
tion"). This equation is well known in the theory of univalent functions (see, e.g., [21]) as 
a differential equation for conformal maps of slit domains parameterized by a parameter 
Q- 

In the calculation below, we closely follow [12]. Let us plug the ansatz (15^ into the 
equation 0131) . Using the chain rule of differentiation, we have: 



ei{w{zi) - w{z2)) 



D{zi)q ■ dq log{e''w{z2)) = - , ^777^^ dt^q 
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Now, using ffTTl) and assuming that 7^ 0, we get the relation 

2 w[zi) - W[Z2) 



which means that the combination 



1 + dqlog W{z) 



viQ) = -wiz)- — — (36) 

I — Oq log W[Z) 

does not depend on z. This imphes the radial Lowner equation fl5^ or, for the Lax 
function z{w, q), 

dzjw) ^ _ ^ viq)+w dz{w) 
d q r]{q) — w dw 

Note that the functions l/w{z) and z{w~^) obey the same Lowner equations (IMl) and 
(1371) respectively. 

Equation (1141) implies that \ri{q)\ = 1. Indeed, with the reduction imposed it becomes 

lo, rq-r-wfi^ r q r w dg\og{w{zi)w{z2)) 
-d,logieM^,))dMeMz2)) = 

or, after rearranging, 

_ 1 + dq\ogw{zi) 1 + dg\ogw{z2) ^ 



l-dq\ogw{Zi) l-dq\0gw{z2) 



that just means that ri{q)ri{q) = 1. In fact this constraint follows already from eq. (136|) : 
let z belong to the image of the unit circle under the map z{w) (i.e., to the Lax contour), 
then complex conjugation of (136|) yields ri{q) = rj~^{q). 

Using the Lowner equations for z{w) and z{w~^), it is straightforward to verify that 
{z{w) , z{w~^)} = 0. As it was argued in section 3, this means a t^-iiidependent relation 
R{z, z) = between the Lax functions z and z which defines a curve F in the plane. Let 
us represent it in the form z = Sr{z). The function Sr is called the Schwarz function 
of the curve F [22]. For the degenerate solution of rank 1 under consideration it is an 
integral of motion. It is easy to see that the canonical transformation (fT9!) corresponding 
to this solution can be written in terms of the Schwarz function as follows: 

-z = l/Sr{z), M = -^I^^M (38) 

The reality constraint fl20l) follows from the identity Sr{Sr{z)) = z obeyed by the Schwarz 
function. 

As w sweeps the unit circle, z{w) sweeps an arc of the curve F (back and forth). The 
arc depends on all the times through q. The function z{w) conformally maps the exterior 
of the unit circle onto the complement of the arc. We see that this map does obey the 
Lowner equation as it must. 

The dependence of q on the times tk, 4 is determined by a system of equations 
of hydrodynamic type. They follow from eq. (fTT!) which can be written as D{z)q = 
— |(1 + dq\ogw{z))dtQq. Using the Lowner equation, we obtain: 

D{z)q = , f , , d,,q (39) 
w[z) -r]{q) 
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From ([5]) we conclude that 

with (pk{ui) as in flTTI) and thus the system of equations of hydrodynamic type reads 

dt,q = MviqMoq , k=l,2,... (40) 

Equations containing t^-derivatives are obtained by complex conjugation. 

At last, it should be mentioned that the chordal version of the Lowner equation (see, 
e.g., [13]), emerges, in a similar way, in the context of the dispersionless KP hierarchy 

mm- 



6 The large N Dyson gas representation of solutions 
to the dToda hierarchy 

In this section, we reconstruct the solutions of the dToda hierarchy (both non-degenerate 
and degenerate) using the Dyson gas representation, i.e., eigenvalue versions of matrix 
integrals for certain models of random matrices. 

Consider the following A^-fold integral over the complex plane: 

JV . J'^ m<n j=l 

where dfi is some integration measure and ^ is a parameter. For dfi = e~T^^'-^'^^d'^z the 
integral is equal to the partition function of the model of normal random matrices with 
the potential 27le J2k ^kz'' — U{z, z) written as an integral over eigenvalues. Equivalently, 
it is equal to the partition function of the system of 2D Coulomb charges interacting 
via the logarithmic potential in an external field (the Dyson gas). It appears that both 
Laplacian growth and growth of slit domains can be simulated by the large N limit of 
this integral. 

The basic fact linking the integral fj^T]) to integrable systems is that for any measure 
dfi (including singular measures supported on sets of dimension less than 2), tn, as a 
function of {tk}, {ik}, is a r-function of the 2D Toda hierarchy with a nonzero dispersion 
parameter proportional to h, i.e., it obeys the full set of Hirota bilinear identities for this 
hierarchy [221 E]- In a slightly different form, this statement first appeared in [23], see 
also [25]. The dispersionless version is reproduced in the large limit such that N oo, 
h 0, and to = hN remains finite. Then r^v generates the dispersionless "r-function" 
(or rather "free energy" ) JF via 

J^(to, {tk}, {h}) = lim (h^ log tn) (42) 

It obeys the dispersionless Hirota relations (see below). 

Second order t^-derivatives of T enjoy a nice geometric interpretation through confor- 
mal maps. This goes as follows. As — > cxd, the integral fHT]) is determined by the most 
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favorable configuration of Zj's, i.e., the one at which the integrand has a maximum. Using 
the electrostatic analogy, one can see that this holds when the points z-i (2D Coulomb 
charges) densely fill a bounded domain D in the complex plane. In terms of the mean 
density of the charges {p{z)) = h i^^6'^'^\z — Zk)^ this domain is characterized by the 
condition 

lim {p{z)) > if z G D and lim {p{z)) = otherwise 

For simplicity, we assume that D is connected. In the matrix model interpretation, this 
domain is called the support of eigenvalues. 

Let w{z) be the conformal mapping function from the exterior of the domain D onto 
the exterior of the unit circle normalized as w{z) = z/r + 0{1) at large \z\ with a real r 
called the exterior conformal radius of the domain D. In [261 [IH] it was shown that the 
function w{z) can be expressed through JF in the following different but equivalent ways: 

rw{z) = ze-^'^'^^'o^ (43) 

rw{z) = z-a- D{z)dt,T (44) 

rw-^{z) = D{z)dt,J^ (45) 

where 

and the operator -D(-2) is defined in section 2. The consistency of these relations follows 
from equations of the dToda hierarchy. The 2D dToda hierarchy can be written in a 
generating form as 

C(..)C(..)^ = log (47) 

^1 — ^2 

- D{z,)D{z^)T = log f 1 - -^-^—^] (48) 



together with complex conjugate equations ( cf. eqs. (1131) . (1141) which are to-derivatives 

of iZD, (BHD). 

We emphasize that all the relations given above hold true for any measure dp in 
(14T|) provided the most favorable configuration of Zj's at — oo is well defined. If the 
measure is smooth, say dp = e~Ti'^^^'^^d'^z, then this construction gives non-degenerate 
solutions to the dToda hierarchy discussed in section 3. It is easy to see that U{z, z) 
is just the generating function of the canonical transformation (|2T1) . hence the notation. 
Singular measures dp lead to degenerate solutions. In particular, one may consider the 
measure supported on a curve F, then the integral ( 14TI) becomes one-dimensional (along 
F) in each variable: 

1 N _ 

= 4 / n - ^-r n e^^^^^"^"-''"^^ m m 

• •'^ m<n j=l 

In the large N limit, the support of eigenvalues, D, is then an arc of the curve F (or several 
disconnected arcs). The function w{z) maps the slit domain C \ D onto the exterior of 
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the unit circle. The choice of the measure supported on a curve means a reduction of 
the dToda hierarchy. A famihar example is the dToda chain, where one may take F to 
be either real or imaginary axis. Consider a general (continuous) curve F infinite in both 
directions. It is clear that w{z) and the Lax functions (the functions inverse to w{z) and 
w{z)) depend on the times through two parameters only. One can set them to be, for 
example, the positions of the two ends of the arc D on the curve F. This is a reduction 
of rank 2. The simplest way to obtain a reduction of rank 1 is to take the measure dfi 
supported on a half-infinite curve starting at a point zq such that the arc D always starts 
at zq as the times independently vary in some open set. 
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